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Abstract 

The quality factor and mode volume of a nanocavity play pivotal roles in realizing the strong coupling interaction 
between the nanocavity mode and a quantum dot. We present an extremely simple method to obtain the mode 
volume and investigate the effect of the slab thickness on the quality factor and mode volume of photonic crystal 
slab nanocavities. We reveal that the mode volume is approximatively proportional to the slab thickness. As 
compared with the previous structure finely optimized by introducing displacement of the air holes, via tuning the 
slab thickness, the quality factor can be enhanced by about 22%, and the ratio between the coupling coefficient 
and the nanocavity decay rate can be enhanced by about 13%. This can remarkably enhance the capability of the 
photonic crystal slab nanocavity for realizing the strong coupling interaction. The slab thickness tuning approach is 
feasible and significant for the experimental fabrication of the solid-state nanocavities. 
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Background 

Photonic crystals (PCs) [1-3] are artificial dielectric nano- 
structures with a periodic variation of dielectric function 
in the length scale of optical wavelength, and provide a 
unique way to control the decay kinetic of the quantum 
emitters inside the PCs due to photonic bandgaps and a 
strong inhomogeneity of electromagnetic fields [4-10]. 
Since many sophisticated and mature fabrication technolo- 
gies developed in micro-electronics and opto-electronics 
can be applied to its fabrication, the PC slab, which is a 
thin semiconductor slab with two-dimensional (2D) peri- 
odicity along the slab plane, has been investigated energet- 
ically in depth both theoretically and experimentally 
[11-15]. Owing to the strong vertical optical confinement 
and the 2D photonic bandgap effect, the overall spontan- 
eous emission rate of the quantum emitter inside the PC 
slab decreases substantially [14]. 

By introducing an artificial point defect into the PC 
slab, the PC slab nanocavity [3] can be formed. The 
point defect traps a localized nanocavity mode, which 
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decays in inverse proportion to the quality factor of the 
PC slab nanocavity. The PC slab nanocavity and a single 
two-level quantum dot can realize the strong coupling 
interaction and thus constitute the solid-state strong 
coupling system (SSSCS) [16]. In this SSSCS, there is 
reversible exchange of a single photon between the 
quantum dot and the nanocavity mode before the pho- 
ton leaks out of the nanocavity. The SSSCS realizes 
many fascinating but genuine quantum behaviors in cav- 
ity quantum electrodynamics [17], e.g., vacuum Rabi 
splitting [16,18,19] and lasing under strong coupling 
[20]. The SSSCS not only provides test beds for funda- 
mental quantum physics but also has important applica- 
tions in quantum information processing [21-23]. 

The realization of the strong coupling interaction re- 
lies on the condition that the coupling coefficient be- 
tween the nanocavity mode and the quantum dot 
exceeds the intrinsic decay rate of the nanocavity [17]. 
To fulfill this condition, a great deal of efforts [24-27] 
have been devoted to design the nanocavities with the 
ultrahigh quality factor and ultrasmall mode volume. 

To enhance the quality factor, various types of the PC 
slab nanocavities have been presented. The prominent 
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types of the PC slab nanocavities with ultrahigh quality 
factor include the PC L3 nanocavity [25] and PC 
heterostructure nanocavity [27]. The PC L3 nanocavity 
is formed by missing three air holes in a line and dis- 
placing several pairs of air holes at both edges of the 
nanocavity, which can increase the quality factor sub- 
stantially by following the principle that light should be 
confined gently in order to be confined strongly [25,26]. 
The PC heterostructure nanocavity is formed by adjus- 
ting the lattice constant of several rows of air holes and 
introducing mode gap difference in the PC slab wave- 
guide, which can obtain unprecedentedly ultrahigh 
quality factor by following the same principle [27]. Obvi- 
ously, this principle requires the elaborately designed 
and optimized PC slab nanocavity with highly fine tun- 
ing of the positions and radii of the air holes around the 
nanocavity center, commonly up to the nanometer scale 
accuracy, which is a great challenge due to the technical 
limits of the semiconductor process [28]. However, the 
effect of the PC slab thickness on the quality factor has 
not been reported. 

Besides the quality factor, another important param- 
eter for the realization of the strong coupling interaction 
is the mode volume of the nanocavity. Traditionally, the 
mode volume is calculated by simulating and then inte- 
grating the electric field distribution of the nanocavity 
mode around the whole nanocavity region [24-26,29] 
(see Equation 6). This is a rather time-consuming and 
difficult task. Obviously, a simple and efficient numerical 
method for the calculation of mode volume is desirable 
and remains a challenge so far. 

In this paper, we present an extremely simple 
method to determine the volume of a nanocavity 
mode and investigate the effect of the slab thickness 
on the quality factor and mode volume of the PC slab 
nanocavities based upon projected local density of 
states for photons [30]. It is found that the mode vol- 
ume monotonously expands with the increasing slab 
thickness. As compared with the previous structure 
finely optimized by introducing displacement of the 
air holes, via tuning the slab thickness, the quality fac- 
tor can be enhanced by about 22%, and the ratio be- 
tween the coupling coefficient and the nanocavity 
decay rate can be enhanced by about 13%. Our work 
provides a feasible approach to manipulate the quality 
factor and mode volume in the experiment. This is 
significant for the realization of the strong coupling 
interaction between the PC slab nanocavity and a 
quantum dot, which has important applications in 
quantum information processing [21-23]. 

Methods 

The optical properties of an arbitrary dielectric nano- 
structure can be characterized by the projected local 



density of states (PLDOS) [30], which is defined as fol- 
lows: 



Z|d.E A (r 0 )| 8((o-a>\), 

x 



(i) 



where r 0 is the location; co, the frequency; d, the orienta- 
tion; and E A (r) and cox, the normalized eigen electric field 
and eigen frequency of the Ath eigenmode of the nano- 
structure, respectively. 

In an ideal single-mode nanocavity without loss, the 
PLDOS can be expressed as follows: 



p c (r 0 ,w,d) = |d.E c (ro)| 2 £(w-6) c ), 



(2) 



where E c (r) and co c are the normalized eigen electric field 
and eigen frequency of the nanocavity mode, respectively. 

Considering the loss, the PLDOS of a realistic single- 
mode nanocavity can be generalized to Lorentz function 
[31] as follows: 



|d-E c (r 0 )| 
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(3) 



where k = cj c I Q is the decay rate of the realistic 
nanocavity with loss and Q represents the quality factor. 
Apparently, when k is infinitely small, Equation 3 of the 
loss nanocavity approaches to Equation 2 of the lossless 
nanocavity. 

For a specific location r 0 and orientation d , the 
PLDOS of the loss nanocavity reaches its peak value p cp 
at o) = co r as follows: 



2|d.E c (r 0 )| 



cp 



(4) 



so we can further simplify the PLDOS of the loss 
nanocavity as follows: 



cp 



k/2 
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(5) 



The mode volume of the nanocavity characterizes the 
confinement and localization of the nanocavity mode 
and is defined as follows [26]: 



V ■ 



Je r (r)|E c (r) 


Vr 


max 


«r(r)|E c (r)| 2 ] 



(6) 



where e r (r) is the relative dielectric constant and E c (r) is 
the electric field of the nanocavity mode. The numerator 
is the normalization factor of the nanocavity mode field. 
The calculation of the normalization factor is rather 
difficult and time-consuming. However, since we can 
directly use the normalized nanocavity mode field E c (r) 
adopted in Equations 2 to 4, we do not need to calculate 
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this normalization factor. With the normalized nanocavity 
mode field E c (r), Equation 6 can be simplified as follows: 



max[£ r (r)|E c (r)| 2 ] 

We assume that s r (r) | E c (r) | 2 reaches to its maximum 
at location r 0m and denote the direction of the vector E c 
(r 0m ) at this location as d m . For most of the PC slab 
nanocavities, r 0m and d m are known before the simula- 
tion. For instance, for the PC L3 nanocavity, r 0m is at the 
nanocavity center and d m is perpendicular to the line of 
centers of the three defect air holes, as will be shown in 
Figure lb. 

By substituting Equation 4 with Equation 7, we can 
obtain the following: 



nK Pcpm E r(*Om) ' 

where p cpm = p c (r 0m ,o) c ,d m ^ is the peak value of the 

PLDOS at the location r 0m along the direction d m . 
Therefore, as soon as the PLDOS at the location r 0m 
along the direction d m is calculated by various numerical 
methods, G) c , k, and p cpm can be determined by fitting 
the PLDOS by the Lorentz function of Equation 5. Based 
on them, we can finally obtain the mode volume of the 
nanocavity by Equation 8 and the quality factor of the 
nanocavity by Q = cj c I k. 

Traditionally [24-26,29], the mode volume of the PC 
slab nanocavity is calculated directly by Equation 6. By 
this method, the electric field distribution of the 



nanocavity mode around the whole nanocavity region 
needs to be simulated and then integrated. This is rather 
time-consuming. In contrast, using our method of Equa- 
tion 8, we can calculate the mode volume simply and effi- 
ciently. We just need to calculate the PLDOS at only one 
known location and along one known direction, which 
make the calculation of the mode volume very efficient. 

As mentioned previously, the realization of the strong 
coupling interaction requires that the coupling coefficient 
g exceeds the intrinsic decay rate of the nanocavity mode 
k. Thus, the most important issue for the realization of 
the strong coupling is to increase the ratio of g/i^Q/ 
\fVo) c [16], which can be obtained using our method. 

Results and discussion 

To compare our slab thickness tuning approach with 
previous air hole displacement approach, we investigate 
the PC L3 nanocavity that was finely optimized by the 
air hole displacement approach in [26], as shown in 
Figure la. 

The 2D PC slab is composed of silicon (refractive 
index n = 3.4) with a triangular lattice of air holes. The 
lattice constant is a = 420 nm. The slab thickness is 
d = 0.6a, and the air hole radius is r = 029a. The PC L3 
nanocavity is formed by missing three air holes in a line 
in the center of the PC slab and can be further opti- 
mized by firstly tuning the displacement A of the first 
nearest pair of air holes and then tuning the displace- 
ment B of the second nearest pair of air holes and, 
finally, the displacement C of the third nearest pair of 
air holes, as shown in Figure la. 
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Figure 1 The structure diagram and nanocavity mode of the PC L3 nanocavity. (a) Cross section on the central plane (z = 0 plane) of the 
PC L3 nanocavity. Gray region is the dielectric slab, and white regions are the air holes. A, B, and C denote the displacements of the first, second, 
and third nearest pair of air holes, respectively. The air holes are moved outward along the x direction, denoted by the arrows, (b, c) E y 
component of the electric field E c (r) of the PC L3 nanocavity mode with the air hole displacements A = 0.2a, B = 0.025a, and C = 0.2a (b) on 
z = 0 plane and (c) on y = 0 plane, respectively. The electric field distribution is normalized by the electric field maximum at the center of the 
nanocavity r 0m = (0, 0, 0). The two dotted lines denote the top and bottom surfaces of the slab. 
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The E y component of the electric field E c (r) of the 
nanocavity mode is shown in Figure lb,c, obtained by 
finite-difference time-domain method [32]. This spatial 
distribution is typical among all the PC L3 nanocavities. 
Obviously, most electromagnetic energy of the nano- 
cavity mode is localized in the three missed air holes 
due to the 2D photonic bandgap effect and is also con- 
fined inside the slab by the total internal reflection. The 
E y component reaches its maximum at the nanocavity 
center r 0m = (0, 0, 0). 

First of all, we focus on the cases where the slab thick- 
ness is fixed at d = 0.6a, and the air hole displacements 
A, B, and C are tuned and optimized in turn according 
to [26]. The PLDOS of the non-optimized and the three 
optimized PC L3 nanocavities are calculated, and the re- 
sults are shown in Figure 2a. Obviously, as the PC L3 
nanocavity is further tuned and optimized, we find that 
(a) the resonant frequency slightly shifts to the lower fre- 
quency, and (b) the decay rate of the PC L3 nanocavity, 
i.e., the full-width at half maximum of Lorentz function 
of the PLDOS, is further suppressed, which leads to the 
remarkable increase of quality factor, as shown in 
Figure 2b. 

However, as the three pairs of air holes near the PC L3 
nanocavity center are further moved outward, the 
nanocavity mode is confined inside the nanocavity more 
and more gently [25], as shown in Figure lb. Consequently, 



the mode volume of nanocavity mode becomes large, as 
shown in Figure 2c. The calculated mode volume of the 
optimized PC L3 nanocavity with air hole displacements A 
= 02a, B = 0.025(2, and C = 02a is 0.0754 \im 3 , which 
agrees well with the reported mode volume as 0.074 (im 3 
in [26]. This excellent agreement validates our method of 
Equation 8 for calculating the mode volume. 

Based on the calculated quality factor, resonant fre- 
quency, and mode volume, we can obtain the ratio of 
g/K, which assesses the PC L3 nanocavity for the 
realization of the strong coupling interaction between a 
quantum dot and the nanocavity mode. As the air hole 
displacements A, B, and C are tuned and optimized in 
turn, g/K is also enhanced remarkably, as shown in 
Figure 2d, which is mainly due to the sharply decreased 
decay rate k of the nanocavity. 

Actually, based on the previous optimized PC L3 
nanocavity with air hole displacements A = 02a, B = 
0.025a, and C = 02a, we can further enhance the quality 
factor by optimizing its slab thickness. We calculate the 
PLDOS of the PC L3 nanocavities with different slab 
thicknesses. The results are shown in Figure 3a. As the 
slab thickness increases from d = 0.5a to d = 1.0a, the res- 
onant wavelength of the PC L3 nanocavity also increases, 
and hence, the resonant frequency decreases substantially. 

As shown in Figure 3b, as we tune the slab thickness, 
the quality factor varies remarkably and reaches its 
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Figure 2 The PC L3 nanocavities with the slab thickness d = 0.6a and different air hole displacements. Including no displacement' 
(denoted as No), 'A = 0.2a' (denoted as A), A = 0.2a, B = 0.025a' (denoted as AB), and A = 0.2a, B = 0.025a, C = 0.2a' (denoted as ABC), (a) The 
PLDOS at the center of the PC L3 nanocavities, orientating along the y direction, normalized by the PLDOS in vacuum as oj 2 I 3n 2 c 3 . (b) The 
quality factor, (c) The mode volume, (d) The ratio of g/K. 
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Figure 3 The PC L3 nanocavities with different slab thicknesses. The air hole displacements are A = Ola, B = 0.025a, and C = 0.2a. (a) The 
PLDOS at the center of the PC L3 nanocavities, orientating along the y direction, normalized by the PLDOS in vacuum as u 2 I 3n 2 c 3 . Each 'vertical 
line' is actually a Lorentz function with small full-width at half maximum, (b) The quality factor, (c) The mode volume, (d) The ratio of g//c 



maximum at the slab thickness d = 0.8a. By the slab 
thickness tuning approach, we can further optimize the 
quality factor from Q = 265,985 for d = 0.6a in [26] to Q 
= 325,121 for d = 0.8a, with increase of about 22%. This 
optimized PC L3 nanocavity with higher quality factor is 
desirable and beneficial to the realization of the SSSCS. 

Along the vertical (z) direction perpendicular to the 
slab plane, the electric field of the nanocavity mode is 
mostly confined inside the slab by the total internal re- 
flection, as shown in Figure lc. Thus, when the slab 
thickness increases from d = 0.5a to d = 1.0a } the 
nanocavity mode is confined inside the slab more and 
more loosely, and hence, the mode volume expands al- 
most linearly along with the increasing slab thickness, as 
shown in Figure 3c. 

As we tune the slab thickness, the ratio of gh< varies 
substantially and also reaches its maximum at the slab 
thickness d = 0.7a. The optimized gh< at the slab thick- 
ness d = 0.7a is about 13% higher than that of d = 0.6a 
in [26]. From Figure 3d, we can notice that there is an 
optimization region for the slab thickness from d = 0.7a 
to 0.8a, in which the ratio gh< varies little. This is very 
beneficial for the experimental fabrication of the PC L3 
nanocavity. 

In a word, the nanocavity mode is confined inside the 
PC L3 nanocavity by the 2D photonic bandgap effect 
along the slab plane and also by the total internal 



reflection in the out-of-plane direction. Thus, as we dis- 
place the air holes near the nanocavity center outwards 
or as we increase the slab thickness, the nanocavity 
mode is confined inside the nanocavity more gently and 
loosely. In this case, the mode volume of the nanocavity 
mode expands, and the electric field maximum at the 
nanocavity center decreases, which results in the de- 
crease of the coupling coefficient g between a quantum 
dot and the nanocavity mode. Since the ratio gh< be- 
tween the coupling coefficient and the nanocavity decay 
rate characterizes the capability of the PC L3 nanocavity 
for realizing the strong coupling interaction between a 
quantum dot and the nanocavity mode, we should pay 
more attention to enhance the ratio gh<, instead of only 
pursuing higher quality factor. 

Conclusions 

In summary, we have presented a simple and efficient 
method based upon the projected local density of states 
for photons to obtain the mode volume of a nanocavity. 
The effect of the slab thickness on the quality factor and 
mode volume of photonic crystal slab nanocavities has 
been investigated, which both play pivotal roles in cavity 
quantum electrodynamics. 

We find that the mode volume is approximatively pro- 
portional to the slab thickness. Furthermore, by tuning 
the slab thickness, the quality factor can be increased by 
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about 22%, and the ratio gh< between the coupling coef- 
ficient and the nanocavity decay rate can be enhanced 
by about 13%, as compared with the previous PC L3 
nanocavity that is finely optimized by introducing 
displacement of the air holes at both edges of the 
nanocavity. Based on these results, we can conclude that 
the optimization of the slab thickness can remarkably 
enhance the capability of the PC slab nanocavity for the 
realization of the strong coupling interaction between a 
quantum dot and the nanocavity mode. The slab thick- 
ness tuning approach is feasible and significant for the 
experimental fabrication of the solid-state nanocavities. 
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PC: photonic crystal; PLDOS: projected local density of states; SSSCS: 
solid-state strong coupling system; 2D: two-dimensional. 
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